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Abstract 

We develop a new and efficient method to systematically analyse four dimensional 
effective supergravities which descend from flux compactifications. The issue of finding 
vacua of such systems, both supersymmetric and non-supersymmetric, is mapped into 
a problem in computational algebraic geometry. Using recent developments in com- 
puter algebra, the problem can then be rapidly dealt with in a completely algorithmic 
fashion. Two main results are (1) a procedure for calculating constraints which the 
flux parameters must satisfy in these models if any given type of vacuum is to exist; 
(2) a stepwise process for finding all of the isolated vacua of such systems and their 
physical properties. We illustrate our discussion with several concrete examples, some 
of which have eluded conventional methods so far. 



1 Introduction 

The issue of moduli stabilisation is one of the most pressing in string phenomenology today. 
Recent progress in this field has resulted in a variety of reasonably well-understood, com- 
pletely stable vacua [H El El IH Ej ■ However, these vacua, for the most part, are not physical. 
Two of the greatest problems with these minima from a phenomenological standpoint are 
that they do not spontaneously break supersymmetry and that they give rise to an anti de 
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Sitter external space. Clearly, if we wish to use such vacua as a starting point for building a 
string theoretic description of our world this problem has to be addressed. In the literature, 
this issue is frequently resolved by employing some kind of "raising mechanism," for exam- 
ple, one based on the presence of anti-branes [3110], or on D-terms |71IH]. In the context of 
a well-controlled supergravity descending from a string or M-theory model, there is, how- 
ever, another option. In general, such theories rich in moduli will have vacua which exhibit 
spontaneously broken supersymmetry and which may be de Sitter - even in the absence of 
raising of any sort. 

Finding such vacua is, however, a prohibitively difficult task using conventional meth- 
ods. Generically, a large number of moduli fields are present in four dimensional effective 
descriptions of compactified theories. These describe such features of the internal space as 
its complex and Kahler structure or the form of some vector bundle, to name but a few. 
Therefore, one is confronted with potentials of supergravity theories as complicated functions 
in an overwhelming number of variables. Minimising such an expression can be beyond the 
reach of conventional techniques. 

The purpose of this paper is to present a novel and efficient approach to the systematics 
of finding such flux vacua. In pedagogical detail we provide two basic tools which make the 
search for these extrema relatively easy. The flrst of these is a simple algorithmic process 
for generating constraints on the flux parameters in the superpotential which are necessary 
(and in some cases even sufficient) for the existence of vacua of any given type. The second 
tool we provide is a completely algorithmic way of finding all of the isolated vacua of a given 
system of interest - including non-supersymmetric vacua of the type described above. This 
tool is based upon a method for splitting up systems of polynomial equations into multiple 
systems of simpler such equalities. Thus, we start with a set of equations which describe all 
of the extrema of the potential and break these up into multiple sets of equations, where each 
of these new polynomial systems describes just one of the loci of extrema of the potential 
(say a single isolated vacuum). In the case of isolated vacua these new equations are so much 
simpler than the original expressions that it is found that one can solve them trivially. For 
example, to entice the reader, the following is one of the systems we discuss in later sections 
where we provide concrete examples of our methods: 

K = -41og(-z(f/-f/))-log(-2(Ti-fi)(T2-f2)(T3-f3)), (1) 
W = ^ [4[/(Ti + T2 + T3) + 2T2T3 - T1T3 - T1T2 + 200] . 

This pair of Kahler potential and superpotential has been obtained in the literature by 
compactifying M-theory on a manifold of SU{3) structure ^9J. We call the associated scalar 
potential, as obtained from the usual supergravity formula, V. Solving for the vacua of 
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this model directly by solving the equations dV = is prohibitively difficult, at least as far 
as non-supersymmetric vacua are concerned. Instead, the method described in this paper 
starts by introducing a polynomial ideal (dV), obtained from the (polynomial) numerators 
of the partial derivatives of V. This ideal corresponds to the algebraic variety of extrema 
of V and can be decomposed into so-called primary ideals P{i) by standard algorithms, 
so that (dV) = Pi n . . . n Pn- Each of these primary ideals corresponds to an irreducible 
variety, or, in physical terms, a single branch of the vacuum space. Indeed, each P{i) is much 
simpler than the original one and can be analysed explicitly in many cases. In particular, 
the zero- dimensional primary ideals which correspond to isolated extrema can be studied in 
detail using methods of real algebraic geometry. Applying primary decomposition to (dV) as 
obtained from the above model (subject to the additional, simplifying constraint Re(f/) = 0) 
leads to the following two zero- dimensional primary ideals: 

{3a;2 = 100, ti = 2x, = a;, ts = x, n = 0, ra = 0, rg = 0, y = 0} , (2) 
{9x^ = 500, 5ti = 2x, t2 = x,ts = x, n = 0, = 0, rg = 0, ?/ = 0} . 

Here, we have defined Tj = Tj + itj for j = 1, 2, 3, and U = y + ix. Thus, by breaking the 
equations up in this manner using the techniques we will describe, we render the problem of 
finding isolated extrema of the potential, including its stabilised vacua, trivial. Even if cases 
were to exist where the simplification were not so drastic, this still would not constitute an 
obstacle for us. This is because we provide, in addition, practical algorithmic methods which 
can extract all of the properties of the vacua from these equations, without ever having to 
solve them explicitly. 

In short, the methods we provide are practical and powerful and make short work of 
finding non-supersymmetric vacua and their properties in these flux systems. In slightly 
more technical language, we propose to re-formulate the necessary calculations arising from 
the extremisation of the potential (and, indeed, extremisation problems at large) in terms 
of algorithmic algebraic geometry and commutative algebra. 

We will show the reader that the flux stabilisation problem generically translates to the 
study of saturation and primary decomposition of certain radical ideals in polynomial rings 
over appropriate ground algebraic fields. This rephrasing is far from a need for sophistry, 
but, rather, instantly allows effective algorithms, most of which have been implemented 
in excellent computer packages such as [Ull E], to be applied. In fact, we show that the 
quantities of physical interest are associated with the real roots of complex algebraic varieties. 
Once the affine variety of interest has been processed using the above complex methods, the 
information of physical relevance can be extracted using real algorithmic algebro-geometric 
techniques. In particular, we make extensive use of real root counting and sign condition 
routines based on the theory of Sturm queries. These algorithms then provide us with 
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tremendous amounts of physical information about the vacua of = 1 moduh theories. 

The methods we present find their most natural application within the context of per- 
turbative stabilisation mechanisms, viz., potentials descending from form fluxes, torsion and 
non-geometric effects. In the interests of brevity and clarity we therefore concentrate on such 
cases in this paper. Practically, if one is interested in completely stabilised geometric vacua 
this would imply the consideration of models in type IIA or G2 structure compactifications 
of M-theory. Non-perturbative effects can however be included in this type of analysis and 
we describe how this can be achieved later on in the paper. 

Our approach is very much in the spirit of [T2j where a programme of systematically and 
algorithmically determining the moduli space of A/" = 1 gauge theories, and in particular 
to look for hidden geometric structure in the MSSM, was initiated. Here, we go one step 
further in our computational capability and utilise versatile and productive algorithms in 
both complex and real computational geometry and ideal theory. 

The paper organised as follows. We begin in Section 2 by translating the computation 
of perturbative moduli stabilisation to one of algorithmic algebraic geometry. The problem 
of finding different vacua, SUSY, non-SUSY, Minkowski, AdS, etc., is classified by the type 
of physical questions with which one is faced. We show in pedagogical detail why one is 
led to the study of ideals, their radicals, as well as primary and saturation decompositions. 
Throughout we will focus on the precise algorithms needed for the investigations at hand 
and how they are used in conjunction with one another. At the end of section 2 we recover 
the physical classification presented at the start in a more mathematical context. It arises 
naturally in the process of organising the problem so that it is susceptible to the methods of 
algorithmic algebraic geometry. 

In Section 3, a first example of the utility of the methods we espouse is provided. Using 
a model taken from the literature on non-geometric compactifications 13 , we show how 
the concepts of resultants and their multi-variate generalisation, as well as elimination-order 
Grobner bases, provide us with various constraints which flux parameters must satisfy in 
such models for there to be vacua with various properties. 

In Section 4, we illustrate the various methods described in Section 2 for algorithmically 
finding flux vacua and their properties. This is achieved by applying our methods to a sam- 
ple of problems drawn from the literature, ranging from compactifications of M-theory to 
type II and heterotic string theories. It is demonstrated that indeed the algorithmic meth- 
ods described constitute a conducive path for research in the field, of diverse applicability. 
Finally, we conclude in Section 5. To make the paper self-contained we have included an 
extensive Appendix as a quick guide, first to algebraic geometry and theory of polynomial 
ideals, and second to the actual algorithms in complex and real geometry and commutative 
algebra used throughout the paper. 



4 



2 Flux Vacua and Algebraic Geometry 



We wish to study four dimensional supergravity theories. In the context of moduh sta- 
bihsation, where the chiral superfields of interest are neutral under any gauge group, such 
theories are specified^ by a Kahler potential and a superpotential W . The K and W 
which arise in such string and M-theory phenomenological contexts are not arbitrary. Both 
quantities generically take on certain general forms which are common to all of the pertur- 
bative stabilisation mechanisms currently being investigated in the literature. As such we 
shall concentrate on theories with this structure. 

First, we require that the Kahler potential be taken as a sum of logarithms of (non- 
holomorphic) polynomials in the fields. This class of theories includes the standard form seen 
in the large volume and complex structure limits of string and M-theory compactifications 
of phenomenological interest. These limits are normally considered in discussions of moduli 
stabilisation so that the use of an effective supergravity is justified, and so that explicit 
polynomial formulas can be obtained respectively. We shall briefly describe how to extend 
our methods to other regions of complex structure space later. A typical form for the Kahler 
potential of such a system is as follows: 

K = -\og{S + S)-\og{d,,k{T+T){T^ + f^){T^ + T^)) (3) 
-log(4fc(Z* + Z%Z^ + P){Z^ + Z^)) . 

Here d and d are constants, which could be related to the intersection numbers of the Calabi- 
Yau threefold and its mirror in the case of an SU{?)) structure compactification without 
intrinsic torsion for example. For the discussion at hand such constants will be regarded as 
mere constant parameters; their origin will not be important. 

Next, we must specify the superpotential W . In the same limits of large complex struc- 
ture, volume and weak coupling we again see a common form arising for the perturbative 
superpotentials which are found in moduli stabilisation contexts. The superpotential takes 
the form of a holomorphic polynomial in the fields. This kind of superpotential includes 
all of the perturbative stabilisation mechanisms known to date: flux, geometrical intrinsic 
torsion, and non-geometric elements in the compactification manifold. For example, the su- 
perpotential obtained for the heterotic string with fluxes on an generalised half-flat manifold 
is given as follows jTH HH] : 

W = -^{eo-^Tpo^) + {ea-^Tpa^)Z^+'-{fi''-^rqt)dabcZ'Z'' 

+ ^(/-«T\°)J,,,Z«Z^Z^ (4) 
o 

^For the reader interested in charged fields we note that D-terms can be included trivially in the discussion 
that follows. For the sake of brevity we shall not, therefore, mention them further. 
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Here the e's and /I's are parameters describing the fluxes present in the compactified space, 
while the p's and g's describe the intrinsic torsion. 

Non-perturbative contributions to the superpotential of course will not take the form of 
a polynomial such as (j^. The simplest implementation of the techniques we will shortly 
describe requires the superpotential to be polynomial in the fields. Given the possibility 
of complete perturbative stabilisation in some models we shall adhere to this case for the 
present. Later, we shall return to the issue of non-perturbative contributions to the super- 
potential where we shall describe how these may be accommodated within the structure we 
advocate. 

Given the above Kahler and superpotentials one can proceed, for uncharged moduli fields, 
to construct the scalar potential from the usual formulas ^H]- The scalar potential is given 
by: 

(5) 

As usual the Da represents the Kahler derivative Oa + Qa^K) and /C^^ is the inverse of the 
field space metric 

^AB = dAdsK . (6) 

Given the above-mentioned forms of the Kahler potential and the superpotential, the po- 
tential is a quotient of polynomials in the fields. This feature, together with the polynomial 
nature of W , will be crucial to the methods which we will utilise throughout this paper. We 
note that the potential can still be written as such a quotient even when raising terms such 
as those added in |lj are included. 

In the problem of moduli stabilisation, we are interested in finding the extrema, and 
in particular the minima, of the potential (0). In addition to the supersymmetric minima 
commonly discussed in the literature, for which DW = 0, this will in general include non- 
supersymmetric vacua. These vacua can be de Sitter or Minkowski even in the absence of 
D-terms or any other "raising" mechanisms. Non-supersymmetric minima of this type are 
not normally considered in the literature as even in simple models they are extremely difficult 
to find - a point to which we shall return shortly. The other extrema of the potential are also 
of some interest. The position of maxima neighbouring stabilised vacua, for example, might 
tell us about which set of cosmological initial conditions will allow the system to obtain the 
stabilised configuration. Likewise, such information can make it possible to estimate the rate 
of decay of a metastable vacuum due to tunneling. 



V 



K 



JC^^DaWDbW -3\W 
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2.1 Classification of the Problem 



For clarity, it is expedient to classify the problem at hand into the following subtypes, each 
of which shall be addressed in turn in the ensuing sections. Let there be n fields indexed by 
i, then, the extremisation problem requires that 

diV = fori = . (7) 

We can classify the solutions to ((Zj) by the amount of supersymmetry they preserve, the 
value of the bare cosmological constant they dictate and so forth. We find it useful to define 
the following four subtypes: 



SUSY, Minkowski 


DiW = 0, Vi, W = Q 


SUSY, AdS 


DiW = 0, Vz, W ^0 


NON-SUSY, Partially F-flat 


DiW = 0, i = 1, . . . ,m < n 


NON-SUSY, Non F-flat 


DiW ^oyi 



Now, recall that our potential is a rational function in the fields. As such, the first 
derivatives of the potential can also be written as quotients of polynomials with a related 
denominator. Physically, we are not interested in the solutions to the resulting equations 
which are given by taking the denominator to infinity. These correspond to the infinite field 
runaways common to these models. Therefore, it suffices to confine ourselves to the cases 
where the numerators of the first derivatives of the potential vanish. 

In conclusion then, all the four subtypes of problems in (jH} deal with the vanishing of 
systems of multi-variate (non-holomorphic) polynomial equations. To further simplify we 
circumvent the issue of the presence of both holomorphic and anti-holomorphic terms by 
substituting the expressions for the fields in terms of their real and imaginary parts. This 
then reduces the problem to that of finding the real roots of systems of complex polynomials. 
It should be noted that the problem can also be reduced to such a form in the presence of 
matter, where one would expand the potential up to some given order in these extra fields 
as usual. 

2.2 Mapping the Problem to Algebraic Geometry 

One can try and analytically solve the equations prescribed in ((7j) and (jH)). This can be 
quickly seen to be impossible in all but the most trivial cases. The reason for this is that, 
even if one of the polynomials is of a sufficiently low degree in a given variable to allow for 
an analytical solution, when one substitutes this solution back into the equations to obtain 
a system for the remaining variables the degree of this system with respect to the other 
degrees of freedom is increased. In a very small number of steps the remaining variables all 
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appear with degree five or higher and the system can not be solved. Indeed, solving systems 
of multivariate polynomial equations is notoriously difficult. 

Numerical techniques do not seem to fair any better. Locating the desired minima with 
such methods is intrinsically difficult due to the shallow nature of the minima and the strong 
features generically present elsewhere in the potential. Furthermore, minima of the types 
desired will in general only appear for certain parameter values and this would result in a 
very laborious system of trial and error attempts to find suitable values. We are compelled, 
therefore, to seek more effective methods. 

Our extrema are defined by the vanishing of a set of complex polynomials in the (real) 
fields. Let us temporarily allow the real fields to take complex values. This results in the 
submanifold of (complexified) field space which corresponds to the extrema being defined as 
the locus where a collection of holomorphic polynomials vanish. This is the definition of a 
complex algebraic variety. The reader unfamiliar with algebraic geometry is directed to the 
Appendices where, to make the paper as self-contained as possible, the necessary concepts 
and constructions are provided. Our moduli stabilisation problem is then to find the loci of 
real roots of a complex variety. As described in appendix A any given affine variety can be 
described by ideals in a complex polynomial ring. The extremisation problem of ((Tj) dictates 
that our variety must be defined by an ideal which is generated by the numerators of the 
first derivatives of the potential V. We shall denote this ideal by (dV). 

As a technical point, multiple ideals describe the same variety. For example, as far as 
the physics is concerned, (x) and (x^) describe the same variety, even though the ideals 
themselves as sets of polynomials differ. To neglect such subtle scheme-theoretic differences, 
one can use the so-called radical ideal, which essentially removes trivial powers of the elements 
of the ideal. We denote the radical ideal obtained from (dV) as a/ {dV). To obtain the latter 
from the former, one can use a standard algorithm J7j as implemented in [TIH ITT] . 

Now that we have stated our problem in terms of algebro-geometrical language we may 
proceed to use some of the powerful techniques which have been developed in that field to 
advance our analysis. For clarity of notation let us first tabulate the key symbols which will 
be used throughout; these will be explained in detail in Appendix A. 

2.2.1 Nomenclature 

• I := (/i, . . . , /„) denotes an ideal generated by polynomials fi, . . . , fn- 

• L{I) denotes the variety corresponding to the ideal / and /(M) denotes the ideal 
corresponding to the variety M. There is reverse-inclusion in the sense that L(/U J) = 
L{I) n L{J) and L{I f] J) = L{I) U L{J). 
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• a// denotes the radical of ideal /. Hilbert's Nullstellensatz is the statement on the 
geometry-algebra correspondence: I{L{J)) = ^/j. 

• The quotient of ideal / by J is denoted {I : J). Closely related is the saturation 

of / by J, denoted as (/ : J°°), corresponding geometrically to the sublocus of L(/) 
which does not intersect L{J). 



2.3 Techniques from Complex Algebraic Geometry 

We have now a mathematical object defining the space of extrema of the potential: it is 
the variety L(a/ {dV)) corresponding to the ideal a/ {dV). This variety is not in general 
irreducible. Physically, this simply corresponds to the fact that the extrema of the potential 
may not be connected into one piece. There may be isolated minima and maxima, loci of 
minima with fiat directions and so on. Mathematically, this means that VW) 

is not a 

prime ideal, but rather collectively contains information about all of the different extremal 
loci, the union of which is the extremal variety. Clearly it would be useful to be able to 
separate out the information about, say, lines of maxima, from that of isolated minima. 
Fortunately, a procedure exists in algorithmic algebraic geometry which does precisely this. 

2.3.1 Primary Decomposition 

It is a theorem that any radical ideal such as y^JdV), as we are working over a polynomial 
ring over the complex numbers, is uniquely expressible as an irredundant finite intersection of 
prime ideals. Each prime ideal corresponds to an irreducible variety and physically represents 
a disconnected locus of extrema. The process of finding these prime ideals is a heavily studied 
subject in algorithmic algebraic geometry and is called primary decomposition. A number 
of algorithms have been developed to perform primary decomposition El 120] • We shall 
make extensive use of the Gianni- Trager-Zacharias (GTZ) algorithm ^HI later on in this 
paper when we come to analyse examples and as such a brief introduction to this is included 
in appendix B. This algorithm has been implemented in fTII by GTZ and Pfister. 

If we denote the prime ideal describing the i-th locus by P{i) then, we have the following. 

^/{^v) = P{i)nP{2)n...nP{k) . (9) 

Here k is the number of irreducible components of the extremal variety - the number of 
different loci. The prime ideals P{i) are in general much simpler objects than the reducible 
a/ (dV). As such this process, even on its own, can be of considerable use in attacking 
problems of our kind. This will be seen explicitly once we move on to describe specific 
examples. 
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In summary, we can split up the extremisation problem of ((Tj) by performing a primary 
decomposition of the radical \/JdV) of the ideal {dV). The four subtypes of the problem 
according to (jH)) can, of course, be treated in the same way and we will shortly demonstrate 
this concretely. 

2.3.2 Dimension and Flat-Directions 

Once we have this series of prime ideals describing the various extremal loci for the potential 
of our flux system we can proceed to extract information about the various extrema. The 
extremal manifold L{^y (dV)), using the reverse-inclusion mentioned in Subsection 12.2. II and 
splits up into unions of irreducible pieces: 

H^/W)) = L{P{1)) U L{P{2)) U . . . U L{P{k)) . (10) 

One of the most important things to know about a given locus of extrema is its dimension. 
Our chief interest will be in minima which are isolated in field space; these are fully stabilised 
vacua. 

For an extremum i to be isolated, the dimension of the corresponding prime ideal P{i), 
(or equivalently the dimension of L{P{i))) must be zero^. Physically, the piece L{P{i)) of 
the vacuum would then consists only of discrete points. In general, the i-th extremal locus 
L{P{i)) will not be zero-dimensional, and will exhibit flat-directions, the number of these 
are obviously dictated by the dimension of P{i). We conclude that for all 

Number of Flat directions of locus i = dim(P(z)) . (11) 

Algorithms have been widely developed for computing the dimensions of ideals. A method 
for testing whether an ideal is zero dimensional, for example, is described in Appendix B. 

Once we know the dimensions of the k prime ideals in the decomposition (jH)) we have 
then obtained significant physical information about our system. For example, if we were 
to find that none of the prime ideals are zero dimensional then that flux system would 
have no completely stabilised vacua without flat directions, either supersymmetric or non- 
supersymmetric. If some of the prime ideals are indeed zero dimensional, and if we are 
only interested in isolated vacua, we can then conflne our attention to this subset of the full 
expansion (jUj). 

^The alert reader may be concerned that we are talking about the dimension of a complex variety when 
physically we are interested in the dimension of the space of real roots. For a real root to be isolated it 
is a prerequisite that the complex dimension of the associated P(i) is zero. If this is not the case we may 
simply vary the real part of one of the unconstrained complex fields. It may be the case however that a zero 
dimensional complex variety has no real roots. This is a question to which we will shortly return. 
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Now, we wish to go on to answer more detailed questions. In particular, we are interested 
in the following inquiries. If an ideal is zero dimensional do any of the corresponding extrema 
correspond to real field values? Are the resulting isolated extrema maxima, minima or saddle 
points? Are the extrema in a well controlled part of field space where we can trust the various 
approximations made in obtaining the low energy effective theory we have been studying? Do 
the extrema correspond to dc sitter, anti de Sitter or Minkowski four dimensional universes? 
Are the extrema supcrsymmetric? To answer these questions we need to turn to the subject 
of real, as opposed to complex, algorithmic algebraic geometry. This is the subject of the 
next subsection. 

2.4 Techniques from Real Algebraic Geometry 

We now have some zero dimensional ideals P{i) at hand. As discussed above, we ultimately 
wish to study the real roots of our polynomial system. We now show that it is possible to 
extract the physically relevant information about the extrema of the potential without ever 
finding the explicit location of these roots, in which we mostly have no interest in any event. 
This situation could be compared to the use of algebraic geometry in describing smooth 
Calabi-Yau compactifications. There, we do not know any explicit metric on the internal 
space yet we can still extract much of the physically relevant information. 

As a brief remark, if we primary decompose over the complex numbers it is always possible 
to trivially solve any resulting zero dimensional prime ideals explicitly for the relevant roots. 
The algebro-algorithmic methods described below are still vital, however, for two reasons. 
First, actual implementations of primary decomposition algorithms normally work over the 
rationals where it is not so clear that finding explicit solutions of zero dimensional primes is 
always possible (although we have found in practice it is for these systems - an unexpected 
bonus!). Second, these algorithms can reduce the number of costly primary decomposition 
calculations wc have to perform in analysing a system. These comments will be illustrated 
concretely in later sections. 

Indeed, each polynomial system P{i), can be, by expanding all of the coefficients into 
their real and imaginary parts (or by working over the rationals from the start - which is 
what we do in practice), turned into a system in . . . ,Xn]- We are thus entering the 
realm of real algebraic geometry. In particular, we need to know about the real roots of 
real polynomial ideals. Much less is known about this field than about its complex cousin. 
However, it turns out that some of the few algorithms currently available furnish us with 
exactly the tools we require to extract what we wish to know. 
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2.4.1 Sign Conditions and Real Root Finding 

We will make extensive use of two kinds of algorithms j25- The first kind allows us to 
compute the number of real roots of a zero dimensional ideal (i.e., it allows us to find the 
number of physical isolated extrema of our potential). The second allows us to compute the 
signs of any given set of polynomial functions on each of the real roots of the system |2Ij . 
by means of a so-called Sturm query. A brief description of how these algorithms work is 
provided in Appendix C. Both of these kinds of algorithm have been implemented in ^T] by 
Tobis 122. 

We proceed then by using the first of these algorithms to find the number of real isolated 
extrema of our potential. We then go on to use the second to extract the relevant physical 
information about these extrema. 

Stability of the vacua: The double derivatives of the potential with respect to the fields 
for the system specified in equations (jH)) and Q take the form of quotients of polynomials 
which make up the Hessian matrix ^ ^i^'' ■ In order to check the character of the extremum 

^ axi axj 

one can compute the characteristic polynomial of this Hessian matrix (which is, in fact, a 
rational function) and focus on its numerator polynomial. We can then form the ideal gen- 
erated by the characteristic polynomial and the zero- dimensional primary ideal, describing a 
solution branch and perform an appropriate series of Sturm queries on its roots. This allows 
you to decide algorithmically whether the extremum is a minimum, maximum or saddle 
point. 

Due to the effect pointed out by Breitenlohner and Freedman [2S] it is necessary to 
determine whether these extrema are de Sitter, anti de Sitter or Minkowski before we can 
say whether they correspond to stable vacua. If an extremum is a minimum or saddle point 
with negative cosmological constant it could still be stable. To discover whether this occurs 
in any given case one must check the sign of a certain set of functions :2^. In fact, as phrased 
in [ini, the bound one needs to test, at the critical point Xq of the potential V, is determined 
by the matrix 



If the eigenvalues of this matrix are all non-negative, then the AdS minimum is stable. 
Indeed, for Minkowski or dS, the positive-definiteness of the Hessian matrix ^J^q^^. suffices 
for stability of the minimum. In our case these tests again all turn out to be quotients 
of polynomials and so this can be achieved with the aforementioned algorithms. We can 
therefore determine how many completely stabilised vacua the system has. 




(12) 
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Validity of the effective theory: For these vacua to be in a regime in which our super- 
gravity description is valid we need the values of certain fields, the size of the internal space 
for example, to be much bigger than 1 - let us say greater than 10. By checking the sign of 
the polynomial t — 10, where t is the field under consideration, we can check whether this is 
the case for each of our stabilised vacua. 

Geometry of the vacua: The potential of the system is, as we have already pointed out, 
a quotient of polynomials. As such to deduce whether our extrema correspond to Minkowski, 
anti de Sitter or de Sitter spacetimes it suffices to again find the sign of the numerator and 
denominator. 

Supersymmetry of the vacua: Another important piece of information to have is whether 
the vacua are supersymmetric or not. The F-terms of our system, given Q and are 
again rational functions and so we can check their sign on each of our stabilised, controlled 
extrema. In particular, the algorithms described in Appendix C will tell us if these polyno- 
mials vanish. We can thus determine which of the stabilised vacua are supersymmetric and 
which are not. 

In conclusion, we can learn essentially all of the important information we require about 
the vacua, both supersymmetric and non-supersymmetric, completely algorithmically, with- 
out ever having to explicitly solve the system. Many of the interesting properties of the 
particle physics associated with each vacuum can also be ascertained in this manner. The 
perturbative contributions to the masses and Yukawa couplings in these models, for example, 
are rational functions of the moduli (in appropriate limits). These points raised above clearly 
constitute a very interesting set of questions. We will now pause our general discussions and 
proceed to show how such questions may be attacked concretely. 
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2.5 Saturations and Classification Revisited 

Having tantalised the readers, we now point out a caveat emptur lest they are overwhelmed 
with optimism. In practice, the above discussion has limits when pursued using the prepack- 
aged implementations of the algorithms available in such computer programs as JOl [TT| . 
Indeed, naive applications of the programs often cause them to struggle, halt, or run out of 
memory. 

There are, luckily, various known tricks for avoiding this set of affairs [21]. These tricks 
all fall under the philosophy of splitting principles and are concerned with splitting the 
problem up into more manageable pieces, even before passing the problem to a primary 
decomposition algorithm. 

One key notion in these so-called splitting principles is the idea of a saturation de- 
composition. In this subsection, we will see how this seemingly esoteric technique precisely 
adapts itself to our goal. A more detailed definition and discussion of saturations can be 
found in appendix A. Briefly, given an ideal / and a polynomial /, the saturation, denoted 
(/ : /°°), is equal to 

oo 

sat{I,f):={I:r)=[j{I:n, (13) 

n=l 

where each (J : /"^) is the quotient of / by f^, which is discussed in detail in Appendix A. 
The point is that the saturation (/ : f°°) corresponds geometrically to the space of all zeros 
of the ideal / for which the polynomial / does not vanish ^. 

We now follow the idea in [23] to utilise the splitting principle. Suppose, for some integer 
/, the following identity holds: 

(/ : r ) = (/ : /) • (14) 

In other words, at some finite / the quotient has removed all powers of / from /. Then, we 
have the following decomposition of the ideal J: 

I = {I:r)n{I,f) , (15) 

where (J, /') is the ideal generated by / together with If we take the radical to neglect 
powers, then we have 

Vl^yiTU^ n y^iU) ■ (16) 

•^In fact to be precise the saturation defines geometrically the closure of the complement of L{ f) in L{I). 
If / is one dimensional then there may be zero dimensional points in the variety associated to the saturation 
for which / = for example. In the bulk of this paper, when we will be interested in using saturations, our 
primary concern will be with zero dimensional ideals where this subtlety does not arise. 
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Geometrically, (fT^ is the split we desire: it says that L(J) is the union of a subvariety 
L(a/(/ : /°°)) where / does not vanish, with a subvariety (/, /)) where / does vanish. 

We pause to ask, what is a good choice of polynomial /, or, iteratively, a set of such /'s? 
In general, finding a non-trivial zero divisor, an element / for which (/ : /) 7^ J, can be 
very difficult. For the problem at hand, however, our supersymmetric theories automatically 
provide the perfect choice! These f 's are simply the F-flatness conditions. Recall that one of 
our problems from (jH)), the partial F-fiat case (which computationally is the most illustrative 
case), is to find the solutions to {dV) such that fi = DiW (or, strictly, the polynomial 
numerators of DiW) vanishes only for a subset of fields i = 1, . . . ,m < n. We therefore, 
naturally, choose each F-fiatness equation as an /, iterating from m + 1 to n. Geometrically, 
we can write this saturation decomposition of the vacuum manifold as: 

L{dV) = L{{dV,fuf2,...,fn))U (17) 

U mOV, /l, /2, . . . , f.-U fi+U . . . , /n) : /D) U 

i 

U L{{{{dV, A, /2, . . . , /,_!, /,+!, . . . , /,_!, /,+!, . . . , A) : /D : /f )) U 

L{{{...{dV:fn...--fn-l)--fn))- 

In words, what this decomposition describes is a classification of the different possible vacua 
according to how many of the F-flatness conditions they obey. Thus the flrst term here is 
simply the supersymmetric vacuum space. The second term is the union of all the vacuum 
spaces for which only one of the F-flatness equations is disobeyed, and so on. Once one 
has broken up the problem in this manner one can go on to apply the analysis discussed in 
previous subsections. 

Therefore, this decomposition is physically intuitive, and natural from the point of view 
of the theory of ideals, as well as being practically useful. The classiflcation (jS)) corresponds 
precisely to (fT7|) . The Minkowski vacuum, for example, would be a subset of the flrst 
term, given by L{{dV, fi, . . . , fn, W)), where the superpotential W vanishes in addition to 
all of the F-flatness conditions. Here, a further simpliflcation can be made; indeed, F-flat 
conflgurations are automatically extrema of the potential in supersymmetric systems. Thus, 
the Minkowski vacuum is then L((/i, . . . , /„, W)), 

If we wish to study a given type of vacuum - be it partially F-flat, non- F-flat or com- 
pletely F-flat, all we have to do is to perform the associated saturation decomposition in 
(I17|) . Working with each of these pieces is much more tractable than working with (dV) 
in its entirety. Indeed, some information can be extracted immediately after forming these 
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saturations. For example, if a given piece in the saturation decomposition has a dimension 
of —1 (this is the convention that the system has no roots) then the associated set of vacua 
are absent in the model under consideration. 

We have come full circle and, in the course of setting up a practical method for finding 
minima, have recovered the physical classification (jH)) in the more mathematical context of 
fITTj) . We shall stop our general discussion here. In the following sections, we will address 
each of the subtypes discussed in (jS}, by illustrating with actual examples taken from string 
and M-theory phenomenology. In these specific examples we will find that our method is 
indeed powerful. Primary decomposition breaks the original extremely complicated sets of 
polynomial equations up into more manageable pieces. The prime ideals containing the 
completely stabilised vacua are so much simpler than the full system that they can be often 
solved explicitly - thus furnishing us with a complete knowledge of the vacua we find. 

Let us then proceed to analyse various parts of this expansion for a variety of models. 
Our aim in doing this will be to illustrate the power of this methodology, as well as to see 
what general statements can be extracted in each case. 



3 The SUSY Minkowski Case and Constraints on Flux 



Let us begin with the case of supersymmetric Minkowski vacua. Here, we are solving for the 
vanishing of the superpotential and its derivatives. In this case, some general theory can be 
developed and general, necessary and sufficient, conditions on the fluxes for the existence 
of such vacua can be derived. Similar constraints can be derived in the other cases but in 
those instances these are only necessary conditions. Necessary conditions for the existence 
of non-supersymmetric Minkowski minima in supergravity have also been given in j25j . 

3.1 Resultants and Diophantine Equations 

Before embarking on a full discussion, let us see what happens if there were only a single 
(complex) field. That is, is a degree n polynomial of a single variable x with integer 
coefficients determined by the values of the fluxes. We are therefore solving the system 




(18) 



W\x) 



tti + 2a2X + . . . + nttnX 



,n-l 



. 
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Already, one can learn quite a lot. We know that two univariate polynomials have common 
zeros iff their resultant vanishes j2E]- Therefore, we require that 

Tes{W{x),W'{x)) = 



( CLn 





det 





V ; 



a„. 



{n - l)a„_i 
nan 



O'n-2 



in 



2)an-2 
l)a„_i 



O'n-2 

{n — 1) times 

{n - 3)a„_3 
(n - 2)a„_2 

(n) times 



ai tto 

a2 ai ao 



2a2 cti 
3a3 2a2 ai 



\ 



= . 

(19) 

In general, the resultant of an order m polynomial with an order n one is homogeneous of 
degree m + n in the coefficients. In other words, for our case, the determinant in ()19p is 
a polynomial in the Oj, of homogeneous degree 2n — 1. This is easy to see. Each element 
in the matrix in ()19j) is either or one of the coefficients. Each term in the determinant, 
when expanded, receives one factor from each column. Any non- vanishing term then has the 
same degree as the diagonal term, which is d'^~^a^, of degree 2n — 1. This seemingly trivial 
observation has interesting consequences. It dictates that the resultant vanishes, if and only 
if the coefficients satisfy a homogeneous Diophantine equation. 

Now, recall that in the general problem of studying the critical points of the (ordinary) 
potential there are holomorphic and anti-holomorphic fields in our defining polynomials and 
we needed to expand them into their real and imaginary components and look for real roots 
corresponding thereto. However, our Minkowski problem is simpler in that we need only 
studying the vanishing of the (holomorphic) superpotential and its derivatives, and it suffices 
to find complex roots of a purely holomorphic polynomial system as above. Therefore, we 
can conclude that a Minkowski vacuum exists iff' the resultant, a homogeneous Diophantine 
equation in the fluxes, vanishes. Of course, there is nothing to guarantee that such vacua 
would be physical in the sense that the values of the real parts of the superfields would be 
large and so forth. To check whether this is the case one would have to utilise the methods 
detailed in the ensuing section. 
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As an illustration, let us present the resultant explicitly for some small values of n: 



n 


resultant 


degree 


1 




1 

X 


2 


0-2 (-a? + 4ao 02) 


3 


3 


^3 (— Qi^ 02^ + 4 oi^ as — 18 qq ai + (4 + 27 03^)) 


5 


4 


a4 (—27 oi"^ 04^ + ai^ (—4 03^ + 18 a2 as 04) — 2 oq ai 03 (—9 02 03^ 
+40 a2^ 04 + 96 ao 04^) + oi^ (02^ 03^ — 4 02^ 04 — 6 oq 03^ 04 
+ 144 ao a2 a^^) + oq (~4 02^ «3^ + 16 02"^ 04 + 144 qq 02 03^ 04 
-128 ao a2^ + ao (-27a3^ + 256 ao 04^))) 


7 



(20) 



It would be interesting to study the solutions to such Diophantine equations. The foun- 
dational work on this subject is laid out in ^57], with some recent surveys and results in 

[211 

3.2 Multi-variate Resultants and an Example 

We have discussed the univariate situation above. What about the general case where there 
is more than one variable? In multivariate examples the equivalents of the resultant of the 
previous subsection can be computed algorithmically using an elimination order Grobner 
basis [nni ■ In other words, there is a systematic method of eliminating variables stepwise 
from an ideal, just like Gaussian elimination for linear systems. A description of the algorithm 
for calculating a Grobner basis in the lexicographic ordering - which is an example of an 
elimination ordering - is provided in Appendix A. This elimination, in the uni-variate case, 
produces the resultant discussed in the previous subsection. 

Algebraically this process takes the intersection / fl C[Xi, X„], of the original ideal 
/ C C[Xi, ...,Xn,ai, ...,am] (where the X's are the variables and the a's the parameters in 
the original problem) with the ring C[Xi, X„] of variables to be eliminated. Geometrically, 
this simply corresponds to the projection of the original ideal on to the subspace of the space 
described by the original ring where the eliminated variables vanish. The resultant conditions 
on the a's are then clearly necessary and sufficient for the existence of a root of / for some 
value of the X's. 

Thus, even in the multivariate case, constraints on the fluxes which are necessary and 
sufficient conditions for the existence of supersymmetric Minkowski vacua can still be found. 
In some cases these can be quite compact in form. In others, however, the resulting constraint 
equations can be quite appreciable in size, as we shall see in a concrete example now. This 
constraint on the practicality of resultants of multivariate systems above a certain level of 
complexity is well known fI7\ . 
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Let us illustrate with a concrete example from the literature. Take equation (2.6) of |13j . 
which presents a non-geometric flux superpotential of the form 



W 



ao — 3air + 3a2r^ — a^T^ 

+S(-6o + 36ir - SfosT^ + 63 r^) 

+3[/(co + (ci + Ci + Ci)r - (c2 + C2 + C2)r^ - Csr^), 



(21) 



with the following constraints on the fluxes. 



ao^s — 30162 + 30261 — 0360 
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(22) 



aoC3 + ai(c2 + C2 



C2) - a2(ci + ci - ci) - a3Co = 

C0C2 ~ c\~\~ CiCi — C2C0 = 

C3C1 - C2 + C2C2 - C1C3 = 

C3C0 — C2C1 + C2C1 — C1C2 = 

C2C1 — C1C2 + C1C2 — C0C3 = . 



C062 - C161 + C161 - C260 
C163 - C262 + C262 - C361 
C063 - C162 + C162 - C261 
C162 - C261 + C261 - C360 



















There are also additional constraints which take the same form as those above but with 
the hats and checks switched around. Various useful pieces of algebraic processing of these 
constraints are provided in ^3] • These relations come from, for example, tadpole cancellation 
conditions and integrability conditions on Bianchi identities. 

Finding Minkowski vacua of this system is then the problem of studying the ideal 
/ = {W,drW,dsW,duW} in the ring C(ao,i,2,3, ^0,1,2,3, £0,1,2,3) [5', T, [/], which is a polyno- 
mial ring in variables S, T and U but with all fluxes treated as parameters (formally, we 
call C(ao,i, 2,3, 60,1,2,3, co,i, 2,3) an algebraic extension of the ground field C). If one uses an 
implementation of the relevant algorithms in a package such as piU ITT] then it is assumed 
that none of the flux parameters vanish. The Grobner basis of / in lexicographic order then 
immediately gives that / has negative dimension. In other words, there are no roots in I. 
This is a quite powerful statement without ever solving for anything, or even imposing the 
constraints ((221): there are no Minkowski vacua for this model, if all of the parameters are 
non- vanishing. 

Of course, some flux parameters can vanish. So let us treat them as variables and place 
I in an elimination order Grobner basis, and eliminate S, T, U to obtain our constraints as 
described above. 

The full result for the superpotential given in ()2ip can be obtained in a matter of seconds^. 
The result is a system of 28 constraint equations which the fluxes must obey. We do not 

"^The best way to achieve this is to homogenise the problem, use a Hilbert driven global elimination order 
Grobner basis calculation, and then dehomogenise again at the end. See |3J for details. 
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present these expressions explicitly here as they amount to 8 pages of expressions in this 
font size comprising of 6 degree 3, 12 degree 4, 8 degree 5, and 2 degree 6 polynomials. 

To provide a concrete result in a presentable fashion, let us simplify by setting, for 
example, all 09,3, 60,3 aiid Co,3 to 1. Indeed, there are still many solutions of (j^^ with 
this choice. Now, treat I as an ideal in C[S', T, f/, ai, 02, 61, 62, ci, ci, ci, £2, C2, 62]. We again 
proceed to eliminate 5, T, U using an implementation of elimination orderings in fOJ . 
We find the following constraints as necessary and sufficient for the existence of Minkowski 
vacua^: 

= 80261 - 3ai62 + a2Ci - 62C1 - aiC2 + &1C2, 

= 276i6ici + %lcl - 27blb2C2 + 362C?C2 - 96?c| - SbiCicl 

-27h\ + 276^ - 276fci - %icl -c\ + 276|c2 + %2cl + c^, 
= 27ai62Ci + — 27aibib2C2 + 9aib2CiC2 — 96162C1C2 + ?>b2c\c2 

-9ai6ic^ - 36iCic| - 27ai6? + 27a2&i - 18ai6iCi - 96? ci - Saicf - 66icf - cf 

+18a262C2 + 96^C2 + 3a2c| + 662C2 + c|, 
= 27aia262Ci + 9a2b2cl — 27a{b2C2 + 9aia2CiC2 — 9aib2CiC2 + ?>a2c\c2 

—9a\c\ — SaiCicl — 27a\bi + 270362 — 9a\ci — 18ai6iCi — 6aic? — 36ic? — 

+9a2C2 + 18026202 + 6a2C2 + 362C2 + C2, 
= 21axa\cx + 9a^c? - 27a?a2C2 + 3a2C?C2 - 9a\c\ - 3aiCic^ - 27af + 27a^ 

-27ajci - 9aicj - cf + 27a^C2 + 9a2c| + 4 ■ 

(23) 

Here, Ci = Ci + Ci + Ci and C2 = C2 + 62 + £2. To see, therefore, whether there are any 
Minkowski vacua for the choice of flux values mentioned above, we need only check whether 
the ideal formed by joining and (j^^ over ground field Z has dimension zero or not. 

In fact, in the system specified in ()2ip and (j^^ the Minkowski vacua always exhibit at 
least one fiat direction even when present. It is easy to show that the curve given below 
defines a fiat direction, in the {S, U) plane, of the potential obtained from (pT|) for any 
Minkowski vacuum. 

- 3ai + 6a2To - Sa-^r^ + 5(36i - 662X0 + 363r2) + 3f/(ci - 2c2ro - Scst^) = (24) 

Here tq is the expectation value of the other modulus in the vacuum. 

We would like to emphasise that constraints on the fiuxes such as those given above can 
be obtained in this manner for any of the cases specified in (jH}. To do this, one simply takes 
the relevant piece in the saturation decomposition (treating parameters as variables) and 
eliminates the fields as above. In other words, elimination orderings can provide us with 

^As before one would have to check whether such vacua correspond to physically acceptable field values 
using techniques presented in the next section. 
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necessary conditions on the fluxes for any type of vacuum to exist. However, in other cases, 
due to our complexification of the real field space in order to make the relevant polynomials 
holomorphic, the resulting constraints are only necessary and not sufficient. This is simply 
because the implied roots of the ideal, if the constraints are satisfied, could correspond to 
complex values for the real and imaginary parts of our complex scalar fields. Such roots do 
not of course correspond to physical vacua. In addition, while supersymmetric Minkowski 
extrema are always minima (with the possibility of fiat directions), other forms of extrema 
can be unstable and therefore not correspond to vacua. 

Having discussed how constraints on fluxes can be derived using elimination orderings 
we shall now resume our main discussion. In the next subsection we revert to the question 
of finding vacua in flux systems according to the methods of Section 2. 

4 Attacking the General Problem and Finding Vacua 

With the above prelude on constraints and the Minkowski case finished, let us systematically 
address the question of finding vacua, including the partially F-flat and Non-F-flat cases, 
using our decomposition methods as discussed in Subsection 12.51 We shall consider the full 
expansion ()17|1 in several examples in this section. 

4.1 An Illustrative Example 

Let us completely analyse a simple first example to illustrate our method. Suppose we 
had a four dimensional = 1 supergravity theory defined by the following Kahler and 
superpotential: 

K = -31og(ri+fi)-31og(r2 + f2), (25) 
W = -Tl - T1T2 - T| + lOTi + lOTa - 100 . 

Even this simple example results in complicated equations. Defining Ti = ti + iri and 
T2 = t2 + iT2, the extrema of the potential V is defined by the following: 

= 25(4 + 4(500 - 280t2 + - IOt^ - 10nT2 - 7tI) + Ati{-UOtl + 74 (26) 
+30(100 + + 2T1T2) + 3t2(200 + + AT1T2 + rl)) - 3(20t^ + 4- 60*2(100 
+2rir2 + ri) + 4{7t^ + 10rir2 + 10(-50 + r|)) + 9(10000 + + 2t^T2 
-lOOr^ + 2rir| + + r^{-100 + 3t^)))) , 

= 25(18rj^ + 27r^T2 + r2(5t^ - 60*2 + ^4 - 12ti(5 + ^2) + Qr^) + ^(-900 + lOt^ 
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+7t^ -6ti(10 + t2) + 27r2)) , 

= -25(3tt + t?(60 - 28^2) -ta- 120^2(100 + 2rir2 + r2^)+t?(-1500 + 560t2 

-37tl + 30t^ + SOriTs + 21 t^) + tl{7T^ + IOT1T2 + 10(-50 + t^)) + 4ti(70t^ 

-45(100 + + 2rir2) - 3t2(200 + + 4rir2 + r|)) 

+27(10000 + + 2tIt2 - lOOr^ + 2rir| + + ri2(-100 + Sr^))) , 

= 25(-60t2(ri + r2) + 5t2(ri + 2r2)+t?(5ri + 7r2)-6ti(2(5 + t2)ri + t2r2) 
+9(rf + 3r2r2 + 3rir| + 2r2(-50 + r|))) . 

Solving this system by conventional means is clearly impossible. According to our discus- 
sions, let us, instead, think of (^^1) as an ideal {dV) G ]R[ti, ^2, ''"i, '^"2]- We perform the 
saturation decomposition of ()17|) and present the components thereof in Table ^ We have 2 
complex F-flatness equations: Fy. = Dt^W = 0, i = 1,2. In the table and the text below we 
expand these into 4 real equations and take Re[Frj] = /i, RefFj-j] = /2, Im[Fj'J = and 

With the table we can begin our analysis. First we break up the ideals listed to extract 
any zero dimensional pieces. This part of the analysis ^ is performed using the factorising 
Grobner basis routine [3T] as implemented in ^T]. Once we have a zero dimensional ideal 
we do not decompose it any further at this stage. Anything which is not zero dimensional, 
however, is primary decomposed to check whether it contains any zero dimensional factors. 
We are thus faced with a list of zero dimensional ideals; on these we check for two conditions 
that they must satisfy if they are to describe physical extrema: 

1. The zero dimensional ideal should have real roots; 

2. The real parts of our original superfields should be greater than 1 when evaluated at 
the extrema. 

These checks are performed using the root counting and sign query algorithms based upon 
Sturm queries as implemented in j23 HI] and outlined in Appendix C. 

The first condition is required because our ring variables correspond physically to the real 
and imaginary parts of the physical fields. The second condition is physically motivated. 
This kind of constraint is enforced in systems descending from flux compactifications so that 
the vacua concerned lie both in the large Kahler and large complex structure limits. Large 
values for the real parts of the equivalent of Kahler moduli in these situations are required 

^This ancillary part of the process is not required in the algorithmisation of the problem of finding flux 
vacua and so was not mentioned in section 2. This is simply a practical point - some initial splitting up of 
the relevant ideals in this manner can make the, already quick, calculations involved much faster. 
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if the effective supergravity descriptions being used in these contexts is to be vahd. Large 
complex structure generically leads to the relevant equations being polynomial in the fields. 
Indeed, more general cases can be dealt with in a similar way to non-perturbative effects 
(whose inclusion will be discussed later), at least in certain limits such as the conifold one. 



Ideal 


Interpretation as Vacua 


Physical ? 


(/l, /2, fs, fi) 


tinnprtiA/m mpf ri c 


Yes 

-L CO 




fk) : fn 


yjGii. ulCLll y ±. LLObyj 


No fi — D or fo — n 


where i j ^ k I 






and k, I 


= 1,.. 








((W/1,/2) 


/3~) 


fn 


partially F-fiat 


No, t2 = 


{{{dVJuf-s) 


fn 


fn 


partially F-fiat 


No, No real roots 


mv,fs,f2) 


fn 


fn 


partially F-fiat 


No, No real roots 


((w/1,/4) 


fsn 


fn 


partially F-fiat 


No, No real roots 


mv,u,f2) 


fn 


fn 


partially F-fiat 


Yes 




fn 


fn 


partially F-flat 


No, ti = 


mdV,f.):f^ 


) : fn ■■ fn 


partially F-flat 


No, No real roots 


where i ^ j ^ k ^ I 






and j, fc, / 


= 1,.. 


■ A 






mdv : /r) : fT) ■■ IT) ■■ fn 


non-SUSY 


No, No real roots 



Table 1: Full saturation decomposition of the vacuum (dV) for the potential V and F- 
fiatness equations fi given in the example (|^. 

Examining Tabled we see that we can confine our attentions to just two terms in the sat- 
uration expansion; the two physical ones, corresponding to the supersymmetric (/i, /a, /s, fi) 
and the partially F-fiat {{{dV, fi, /a) : fn '■ fn extrema; both are AdS. Indeed, as well as 
simplifying the analysis this allows us to make quite general statements. For example, all 
non-supersymmetric vacua of this system are partially F-flat with always being zero in 
the vacuum. We proceed to study the two physical extrema in detail. 

We perform the primary decomposition of (/i, /a, /a, fi) using the algorithm due to GTZ 
as implemented in [TT]. The result contains 6 factors all of which are of dimension zero. 
Of these two have 1 real root, one has 2 real roots and 3 have no real roots. Of the 3 
factors having real roots only the single factor with 2 roots is such that the real parts of the 
superfields are valued larger than 1 in the vacua. Thus the physical supersymmetric vacua 
of the system are given by the roots of the following ideal: 

(ti-5,t2-5,ri-r2,9r2-175) C (/i, /s, /a, /4) • (27) 
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Now, for the partially F-flat component, the primary decomposition of {{{dV, f^, : 
fs") '■ fi°) contains 3 factors all of which are again zero dimensional. Of these factors two 
have 2 real roots and one has no real roots at all. Of the two factors with real roots there is 
only one root in one of the factors for which the real parts of the superfields are both greater 
than 1. This is one of the two real roots of the following polynomials (the one for which t2 
is positive): 

(ti - h, 21tl - 20h - 900, n, T2) C {{{dV, U, h) : f^) ■ /D • (28) 

One can ask more about the properties of the above vacua, again using Sturm queries as 
described in Subsection 12.41 and Appendix C. We find that the non-supersymmetric vacuum 
described above is not a local minimum but a saddle point by testing the signs of the 
second derivatives of the potential. Furthermore, the vacuum in question does not obey the 
Breitenlohner-Freedman bound (fT^ and so this vacuum is not stable. Of course, in this 
case the resulting ideals that need to be considered have been rendered so simple by the 
decomposition process that one can simply find the roots of the polynomials in the prime 
ideals analytically. This is in fact generically the case in these flux vacua systems and is 
simply a consequence of the fact that prime ideals tend to take a simple form. 

Solving ()28|1 to find the position of the non-supersymmetric vacuum we obtain ti = ^2 = 
^(1 + Vl90), Ti = T2 = 0. Plotting the potential about this point we can therefore provide 
a check that our method is functioning correctly, as is shown in Figure ^ 

We have presented this example with three main goals in mind. The first is simply to 
give a clear, simple example of the general discussions given in Subsection El The second 
is to demonstrate that this method is practical and powerful. We reiterate that, in the 
system defined by (j^ . we have found all of the isolated vacua of the system. It turns out 
in this case that there are three - two supersymmetric and one non-supersymmetric. To 
find the non-supersymmetric vacua given above and show that it and the F-fiat solutions 
are the only such extrema present in the system using more conventional methods would be 
prohibitively difficult analytically. One would have to find all of the solutions to a system 
of 4 coupled quartics, even in this simple example. Finally, the third goal is to show that 
non-supersymmetric vacua of such systems do exist, even in the absence of D-terms. 

4.2 Examples from String Constructions 

Having whetted the reader's appetite with our toy example, we shall now delve into some 
systems which have been obtained in the literature in the context of string and M-theory 
compactifications to four dimensions. Despite the complexity of the equations which appear 
in these contexts, large portions of the saturation expansion (and in some cases all of it) can 
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Figure 1: The non-supersymmetric vacuum for the supergravity theory specified in (j^EI) 
for our toy example. The fields are Ti = ti + iri and T2 = ^2 + ^^2- We have plotted 
the potential in two slices through field space, viz., ti-t2 and ti-ri. A shift in the V axis 
of 2.09279725 x 10^^ has been performed so that the very shallow vacuum can be plotted 
effectively. 

still be analysed very quickly indeed. In what follows we shall first give a simple example 
from heterotic string theory. We shall then consider an example from M-theory where all of 
the moduli of the system can be stabilised perturbatively without recourse to non-geometric 
spaces. 

4.2.1 A Heterotic Example 

Let us begin with a heterotic theory compactified on one of the SU{3) structure manifolds 
considered in jTU [T3]. Of course, in a heterotic model the dilaton is unstabilised in the 
absence of non-perturbative effects. We shall therefore just consider the stabilisation of the 
analogues of the Kahler and complex structure moduli. In ignoring the dilaton in this manner 
the only modification to the proceeding formulae is that the — 3|iyp term in equation ^ 
becomes — 2|iyp due to a cancellation with the dilaton's F-term. For the Kahler potential 
and superpotential we have [TU [T3] : 

K = -31n(r + f) -31n(Z + Z) (29) 

W = i{i + ieT) + {t + ipT)Z + U^i + iqT)Z^ + ]:{p + irT)Z\ 

I 

where T is the Kahler modulus and Z, the complex structure and ^,r,e,q, fi,p, p,e are 
parameters characterising the fiux and torsion on the internal space. These parameters 
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satisfy the following constraint: 

— eq + pp — pe = . (30) 

As an example let us make the following parameter choices^: 

^ = -13,r = 0,e=-4,g = 2,/x = 2,p=l,p=5,e = -7. (31) 

This gives rise to the following equations for the extremisation of the potential, where 
we have defined T = t + ir and Z — z + i(: 

= Afz^ - Ut^z^ - 25^6 + GOtz^C - 48t'-z\' - JBz^C + 120tz\^ - At^"" - 36T2ct32) 
-75z^C^ + 60tC^ - 25C^ + 24tz^ + 48t^z^C - 60z^C + 36tz^C^ + 8t^C^ + 72t^C^ 
-120^\^ + 12rC^ - 60C^ - 108t^^^ - Ut^z"^ + 360te^ - 12/ + 144rz\ - 60t^C^ 
-540t\2 - 288^2C^ + 636tC^ - 276C^ - 96tz'^ + SGi^C + 504t\ - 192^2C + 1152rC^ 
-1068C^ - im^ - 1764t2 - 192z^ + 1080tC - 1512C^ + 6552t - 3744C - 6084, 

= 2rz^ - + 8rz\^ - 10z\^ + 6rC^ - 5C^ - 2/ - 8t/C - S^C^ - 12tC^ - 
+2tz^ - l2z^C + 90rC^ - 53C^ + Sz^ - 84rC - 96C^ + 294r - 90C - 546, 

= -4iV + 4TV + 25^6-20T^^C-16TVC' + 25A' + 40T^'C'-12tV 

-36t2C^ - 2bz'^C^ + 60tC^ - 25C^ - 8tz^ + 16t2/C + 2Qz\ + 12tz'^C'^ + 24i^C^ 
+72r2(3 _ 40^2^3 ^ ^2tC - 60C^ - lOSt^/ - 4r2z2 ^ 4/ ^ 48r22C - ISOt^C^ 
-540r\2 - 96^2C^ + 636tC^ - 276C^ - 2,2tz^ + im'^C, + 504r2C - 64^^C 
+1152tC^ - 1068C^ - 588t^ - 1764r2 - 64^^ + 1080tC - 1512C^ + 6552r - 3744C 
-6084, 

= -lOr/ + 16r^z\ + 25/C - 60r2^C^ + 8t\^ + 24x^(3 + 50^2^3 _ ^^^^a ^ 25C^ 

-8t2/ + 10/ - 12t/C - ^2t^C - Ser^C^ + 60/C^ - SrC^ + 50C^ - 24tz^ + mt^C, 
+180t^C + 96^^C - 318tC^ + 184C^ - 28i^ - 84t^ + 32z^ - 384tC + 534C^ - 180t 
+504C + 624 . 

^We could in principle avoid choosing parameters by working over an algebraic extension of the base 
field (essentially allowing polynomials with parameter coefBcients) , as was done in the Minkowski example 
in Section 3. Such a calculation would be expensive however. As such, given that the parameters in such 
models are quantised in any case, it is quicker to scan through a given set of values for the fluxes, and 
to automate the following calculations. The calculations involved here are sufficiently quick that this is a 
practical possibility. 
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The algebraic variety defined by these equations is reducible. First of all, we break up the 
variety according to the saturation expansion. Despite the fact that we are dealing with 4 
coupled sextics in 4 variables we can calculate all but the final term (the completely non-F- 
flat case) in (fTTjl for the saturation decomposition extremely quickly. The final term takes 
longer to complete and so we will omit it in what follows. Having obtained the various terms 
in the saturation expansion we go on to study each in turn. 

We again use a mixture of the factorising Grobner basis routine coupled with the GTZ 
primary decomposition algorithm, as implemented in |TTj, to break up the varieties. To find 
out which of the resulting zero dimensional irreducible ideals admit real roots we use the 
appropriate Sturm query algorithms. We also study various sign conditions evaluated on 
these real roots and only keep those vacua for which Re(T), Re(Z) > 0. The only physical 
vacuum that is present is the supersymmetric vacuum which was found in J3] , there are no 
partially F-flat vacua in this system. As such we shall move on to some more complicated 
cases with the aim of finding some non-supersymmetric extrema. 

4.2.2 An M-Theory Example 

Let us now look at another interesting example taken from M-theory. In particular we would 
like to consider a case where all of the moduli are perturbatively stabilised. We will return 
to the question of non-perturbative contributions to the superpotential in the next section. 

One possibility from the literature would be type IIA string theory compactified on an 
orientifold of the ^^^z^ orbifold in the presence of fluxes and torsion, as described in |32| . 
In particular, in their subsection 5.3, these authors provide a choice of fluxes which results 
in a completely stabilised supersymmetric vacuum. 

If one analyses this system using the methodology we have been describing in this paper 
one instantly finds that (dV) for this system contains no zero dimensional ideals at all in 
its primary decomposition. In other words there are directions in field space for which this 
potential is completely fiat. Once the presence of such a fiat direction has been indicated by 
this formalism it is easy to spot it explicitly in the potential - in this case it corresponds to 
a linear combinations of some of the axions of the theory 

Thus, although there is a stable supersymmetric vacuum (supersymmetric configurations 
of this kind automatically obey the Breitenlohner Freedman bound (fT^ ). this system is 
perhaps not of such strong interest for us. For example, there is no hope of finding stable, 
non-supersymmetric, isolated vacua in this model. As such we shall move on to consider 
another possibility. 

An example of the kind we would like, better suited to our purposes, is furnished by [Hj. 
These authors consider compactifying M-theory on the coset ^m^^jfrrfr- This is a manifold 
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of SU{3) structure. The resulting four dimensional supergravity theory is described by the 
following Kahler and superpotential [Oj: 

K = -41og(-2(f/ - U)) - log(-z(Ti - fi)(T2 - f2)(T3 - fs)) , (33) 

W = ^ [4f/(Ti + T2 + T3) + 2T2T3 - T1T3 - T1T2 + 200] . 
v8 

To give an idea of the complexity involved in a case such as this we note that the potential 
takes the form: 

^ = 256tT7i7i^ (40000 + tlrf - 400rir2 - 4tlTiT2 + Atlrl + rfr^ - 400rir3 + 800x2X3+ 
2x^x2X3 - 4rir|r3 + rfrl - AtiT^tI + 4r|r| - 24t2t3a;^ + ^4%^ - 24ti(t2 + ^3)3^^ 
+4rf + StiT2X^ + 4r|x^ + 8x1X3^^ + 8x2X3^^ + 4r|x^ + 1600riy - 8t^ri?/ 
+1600x2?/ + mlr^y - ^r^r^y - Snr^y + 1600r3y - Sr^r^y + 16r|r3y - 8rir|y 
+16r2r|y + Wtly^ + IGrfy^ + 32TiT2y^ + 16r|?/2 + 32rir3i/2 + 32T2T3y^ + 16r||/2 
+t?(ti + tl + r| + 2r2r3 + r| + 4^^ - 8x2?/ - Sr^y + IGy^) + ^2(4^2 ^ ^2 _ 4^^^^^ ^ 2y) 
+4(r| + x2 + 4r3y + 4|/2)) , 

(34) 

where we have defined the component fields by Tj = —itj + Tj for j = 1, 2, 3, and U = —ix + y. 

To obtain the equations for the extrema of this potential we must then take the derivatives 
of this expression with respect to all 8 fields and set them equal to zero. The result is 
somewhat lengthy and so we shall spare the reader the explicit full set of conditions for the 
extremisation of this potential. To solve these equations using normal techniques we would 
have to solve 8 coupled equations in 8 variables with each equation involving a quotient of 
a fourth order and seventh order polynomial, clearly an impossible task (even for packages 
such as Mathematica or Maple). 

However, using our saturation and primary decomposition techniques, the problem is 
much more tractable. Now, in the interests of showing the diverse manners to which our 
methods can be applied, we will present a slightly different analysis for this system. It may be 
the case that one wishes to examine vacua with certain physical properties besides a specific 
degree of F-fiatness. For example, one can ask if there are any vacua for any particular field 
values; for instance, say y = 0. In terms of the variety being considered this is associated 
with the ideal which is generated by dV and the monomial y. This system, which would 
still be prohibitively difficult to solve with more conventional techniques, is well within the 
capabilities of our algorithmic techniques on a desktop computer. The search for such vacua 
might be physically motivated in many ways. For example, one may wish certain axions 
in certain models to vanish in the vacuum in order to agree with a small theta angle in 
a desired target theory. Since we are using the power of this formalism to look at stable, 
non-supersymmetric vacua, demanding such physical inputs hold true is now a reasonable 
thing to do. 
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Again, using a combination of factorising standard basis, GTZ primary decomposition 
and Sturm query algorithms to decompose and analyse the ideal {dV, y) one obtains a de- 
composition involving 16 factors each of which may be made up of numerous prime factors 
themselves. Many of the resulting prime factors are zero dimensional but only two have real 
roots for which the real parts of all of the superfields take values greater than 1. 

As before, the prime ideals which we have extracted from the overall problem to describe 
these isolated loci are so simple that we can solve them explicitly to find the extrema. 
These turning points are described by the following ideals (the generators of which should 
be compared in complexity with the first derivatives of equation ()34|) ): 

/i := {3x'^ -100,ti-2x,t2-x,ts-x,Ti,T2,T3,y), (35) 

I2 := {dx"^ - 500,5ti-2x,t2- x,ts- x,Ti,T2,Ts,y) . 

The simplicity of these equations shows us how useful this procedure is. In separating out the 
ideals that describe the isolated extrema in which we are interested from all of the rest of the 
turning points we have vastly simplified the discussion of the minima - in this case rendering 
it rather trivial. The physical root of Ji is simply the supersymmetric vacuum of the system. 
This reproduces the result found in The physical root of I2 is an isolated extremum of 
the system which is non-supersymmetric and anti de Sitter. These two constitute all of the 
isolated extrema of this system which obey the physical constraint we have imposed. The 
SUSY extremum is Breitenlohner-Freedman stable while the non-SUSY one is not. 

We see that the plots of Figure |21 confirm all of the features of the non-supersymmetric 
extremum that our algorithmic algebro-geometric procedure rapidly predicted. We have also 
calculated a large part of the saturation expansion (fTTj) for this case. We do not however 
find any interesting extrema beyond those described above and so shall not explicitly present 
this analysis here. 

5 Conclusions and further work 

This paper was concerned with the problem of finding vacua of four dimensional supergravi- 
ties describing flux compactifications. After presenting a natural classification of such vacua 
we have provided two primary results within this context. 

First, we have described a practical, algorithmic method for generating constraints on 
the flux parameters in the superpotentials of such systems. We emphasise again that these 
constraints can be derived as necessary conditions for the existence of any given kind of 
vacuum. In the case of supersymmetric Minkowski vacua this result is even more powerful. 
For these special vacua the constraints we have provided are both necessary and sufficient 
for the existence of such extrema. 
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Figure 2: The non-supersymmetric extremum corresponding to the ideal I2 in (jSSl), for the 
supergravity potential specified in (jnH). The fields are Tj = —iti + Ti and U = —ix + y; 
we have here plotted the slices in {ri^t-^) and {t2,y) coordinates. A shift in the V axis of 
4.07 X 10~^ has been performed so that the very shallow vacuum can be plotted effectively. 

Second, and perhaps more importantly, we have outlined a completely algorithmic method 
to find all of the isolated vacua of such systems, be they supersymmetric or not. In addition 
to the vacua themselves the methods we have described enable us to algorithmically find 
most of the quantities of physical interest associated with them. This includes the degree 
of supersymmetry they preserve, their stability as well as particle physics properties such as 
the Yukawa couplings in the matter sector. 

What we have done is to map the extremisation problem to the language of algorith- 
mic algebraic geometry and, in particular, of ideal theory and commutative algebra. This 
is not simply a hypothetical discussion. Using recent advances in computer algebra, the 
methods we present are powerful and allow us to solve, within seconds on an ordinary desk- 
top computer, problems which are simply impossible with conventional techniques. We have 
demonstrated in concrete examples the efficiency with which our algorithms can find isolated 
non-supersymmetric extrema in actual systems directly derived from string and M-theoretic 
compactification. 

One obvious extension of the work presented here would be the inclusion of non-perturbative 
elements in the superpotentials considered. There are several ways in which one might do 
this. The simplest way to proceed would be to simply introduce extra 'dummy' variables 
to represent any exponential functions that appear. One would then have as the desired 
vacuum space an algebraic variety, as described in the bulk of this paper, intersected with 
an exponential equation - that defining the dummy variable. This would enable one to bring 
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the full power of algebraic geometry to bear on the difficult part of the problem. Another 
possibility would be to fix field values at some desired values and then solve the system to 
see what fiux values are required to give stable vacua. In other words, we can solve for a 
set of the parameters rather than the fields. If this is performed carefully this will result 
in an algebraic variety, with the fiuxes as the variables, as the object to be analysed. This 
approach seems to be more difficult to pursue, however, due to the quantised nature of the 
fiux parameters in these systems. 

In any event, in this paper we have restricted ourselves to perturbative superpotentials 
where the methods we have outlined find their simplest application. Such superpotentials can 
result in stabilisation of all of the moduli in geometric IIA and M-theory compactifications. 
Non-geometric compactifications (which give rise to a perturbative superpotential) can give 
rise to stabilised vacua in the other string theories as well. Perturbative vacua are interesting 
as they are on a somewhat firmer footing than their counterparts which rely on a mixing 
of perturbative and non-perturbative effects. One reason for this is that in such mixed 
scenarios one relies on a play off between the two types of superpotential contribution to 
obtain a vacuum. Although such playoffs are theoretically possible with the rest of the 
infinite series of non-perturbative corrections being negligible, such a situation is dependent 
on, for example, a very large coefficient appearing in front of the exponential terms. There 
is no reason to believe that such a coefficient would arise in any given model. 

As a side comment we note that all of the non-supersymmetric vacua we have found thus 
far in any model have been partially F-fiat. 

Further extensions to this work are clear and numerous. As well as the inclusion of 
non-perturbative effects mentioned above one could consider improving the algorithms used 
and their application to the problem at hand. One possible such direction of improvement 
would be to construct a method for performing the calculations over a finite field and then 
separating the spurious results from the physical ones. Grobner basis calculations over finite 
fields can be much faster than those over the rationals. 

Finally, pushing these methods to their natural conclusion, one could imagine a com- 
pletely automated algorithmic approach to extracting the phenomenological physics from 
four dimensional descriptions of string compactifications. Once a four dimensional effective 
theory is derived we have shown that we can scan the vacua of the system and their prop- 
erties algorithmically - searching for appropriate minima with which to describe our world. 
Due to the complexity of these problems [3^ such a program of research would have to be 
guided by physical insight. 
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APPENDIX 

A Rudiments of Computational Algebraic Geometry 



Our computations throughout this paper have relied heavily upon techniques and algorithms 
in algebraic geometry, which may not be entirely familiar to all researchers in the field. With 
this appendix, which will provide a glossary on the key concepts used, we wish that the 
subsequent self-contained nature of this paper may serve the incipience of such methods into 
the study of flux vacua. Detailed exposition can be found in the texts of |2H1 Ell, whose 
emphasis is on the theory, of j30i|, on the computation, of jTUl CH EI], on the practically 
implemented algorithms, as well as of [12] , on a parallel application to Af = 1 gauge theories. 

Algebraic Varieties and Ideals: The problem of finding the vacua of our concern, as 
we stated earlier, is the problem of finding the set M of simultaneous zeros of a system of 
polynomial equations in variables xi, . . . , Xn- Such a set M is an afRne algebraic variety. 
In the language of commutative algebra, in which actual algorithms are always phrased, 
this set is seen as the loci of roots of an ideal I{M) in the ring R = C[xi,X2, ■ ■ ■ , x„] of 
polynomials in Xi with coefficients in C. 

Briefly, recall that a ring is roughly a set with addition (and its inverse, subtraction) and 
multiplication, but no division. Indeed, the sum, difference and product of two polynomials 
remain a polynomial while the ratio does not. An ideal is a subset, which, when multiplied 
by any element, remain in the subset. To intimate the relation between the algebraic object, 
viz., the ideal / and the geometric object, viz., the variety L, the standard notation is to use 
/(L) and L{I) when they correspond. 

To be explicit, we use {fi,...,fk) to denote the ideal of generated by the polynomials /«, 

i.e.. 



for polynomials h^. In this notation, addition and multiplication between two ideals is easily 
defined as the addition and multiplication of all combinations of the generators. Quotients 
will be defined shortly. 




(36) 
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Radical Ideals: Next, I{M) can contain more information than is physically needed. 
Multiplicities in the roots describe the same set of points. Recall the example in the text: 
X = and = describe the same set of points even though (x) and (a;^) are two different 
ideals. This ambiguity is resolved by defining the radical -\/7 of the ideal of / in a ring R: 

77 := {r e i?|r" G / for some n G Z+} . (37) 

The Hilbert NuUstellensatz then states that, for any ideal J, the ideal I{L{J)) corresponding 
to the variety L(J) whose points are determined by J is equal to the radical ideal In 
other words, the radical ideal is the "minimal" ideal corresponding to the variety M which 
drops all the redundant information on the multiplicities of the zeros. Thus, we can refine to 
the study of the radical ideal ^JI{M) corresponding to our zero-set M. Popular algorithms 
which perform this step can be found in ^Tj and are implemented in [TIH ITT] . 

Primary Decomposition: The radical ideal \J I{M) may still be reducible in the sense 
that the variety xy = 0, for example, clearly consists of two irreducible components x = 
and y = 0. To obtain the elemental constituents of a/ I{M) we must then decompose it 
into prime ideals, ideals p for which (just like a prime number), ab & p implies that a & p 
or b E p. Such a process is called primary decomposition®. The theory was originally due 
to Lasker-Noether, with the first algorithm by Hermann. Today, it constitutes one of the 
most exciting areas of research in computer algebra, with popular algorithms by Shimoyana- 
Yokoyama, Eisenbud-Huneke-Vasconcelos, and Gianni- Trager-Zacharias as implemented in 
Pm ITT] . We shall describe the last of these algorithms, which we have used throughout this 
paper, in some detail in Appendix B. We therefore have the decomposition of ^JI{M) as the 
finite intersection of prime ideals P(i), i.e., ^J(M) = ^iP^i). 

Real Roots: After decomposing into irreducible components, one can then compute the 
dimension (corresponding to the number of fiat directions) of each piece P{i)- A method for 
checking whether an ideal is zero dimensional, for example, is briefly described in appendix 
C. In the case that the ideal P{i) is zero-dimensional, the component corresponds to no more 
than a (discrete) set Si of points. Physically, this means that this component of the vacuum 
has been completely isolated. One could determine the cardinality of Si (the number of roots 
of the polynomial system); this is called the virtual dimension of the zero- dimensional ideal 
P{i). In particular, we are interested in the set of real roots, which is a special subset of Si. 

^Strictly, irreducible varieties correspond to primary ideals which are ideals / for which ab E p implies 
that a G p or 6" e p for some integer n, a weaker condition than primality. However, since radicals of 
primary ideals are prime and we are already starting with a radical ideal, it suffices to study the stronger 
condition of prime decomposition. 
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Algorithms have been developed to deal with real roots [21]. We shall discuss this further in 
appendix C. 

Quotients and Saturations of Ideals: During the course of our analysis we need the 
concepts of saturations and quotients of ideals. An ideal quotient of an ideal I <Z R with 
respect to / G -R is simply defined as follows: 

(/ : /) := {g E R\gf E 1} . (38) 

In general, the quotient (/ : J) of an ideal / by an ideal J is the set of elements g E R 
such that g ■ J is contained in I. The definition of a saturation of an ideal is then a simple 
extension of this idea: 

oo 

(/ : D ■■= {g e R\gf E /, for some N E Z>o} = [j I : P . (39) 

n=l 

The second equality is important and is the origin of the infinity in the notation: saturation 
quotients out all powers of /. Geometrically, this means that L{I : f°°) corresponds to the 
subvariety of L{I) for which / 7^ 0. 

Quotient Rings: The last concept that we shall require, for use in later appendices, is that 
of a quotient ring. For an ideal / in a ring R the quotient ring R/I is simply defined to be 
the set of all elements in R where two elements are regarded as equivalent if their difference 
is an element in /. Physically the quotient ring corresponds to the set of all polynomial 
functions where two functions are only regarded as different when they take different values 
on the locus L{I) which is defined by the ideal /. 

Grobner Basis: The first step in almost all algorithms in computational algebraic ge- 
ometry is to place the generators of the ideal of multi-variate polynomials into a so-called 
Grobner Basis. This is a generalisation of Gaussian elimination for a multivariate linear 
system to general polynomials. 

In computational algebraic geometry, the Grobner basis is determined by (modifications 
and improvements of) Buchberger's algorithm (see for example |SI]). The Buchberger algo- 
rithm proceeds as follows. Start with an ideal /. 

1. Set Q = generators (/). 

2. For any pair of polynomials A, B E Q form the 5* polynomial (described below). 

3. Reduce the S polynomial with respect to Q. 



34 



4. If the reduction is non-zero add the result to Q. 

5. Repeat from step 2 until all pairs of polynomials in Q give S polynomials which reduce 
to zero. Q is then the Grobner basis. 

In the above one needs to understand the process of reduction and what an 5* polynomial is. 
Both of these concepts rely on the introduction of monomial orderings. An ordering > is 
simply a rule which allows us to unambiguously compare any two monomials in the variables 
and say which one is higher in a list of all monomials. For example the Lexicographic ordering 
with respect to the variables a, b, c just says that monomials are ordered, firstly according to 
the power of a they contain (highest first), then according to the power of b and finally that 
of c. So, for example, a'^bc would be ordered higher than ab'^c^. 

The reduction process of polynomial A relative to polynomial C is then simply as follows. 
We subtract some (possibly monomial) multiple of C from A in such a manner as to cancel 
A's leading term with respect to the ordering >. If the leading term can not be canceled in 
this way A is simply left alone. 

The S polynomial of two polynomials A and B is simply given as follows. Multiply A 
and B by the lowest degree monomials possible so that the leading terms of the two results, 
A' and B', become equal. One then simply subtracts one from the other, so that the leading 
terms cancel: S = A' — B' . 

Grobner bases have many uses, some of which we shall encounter later in these appendices. 
One particularly useful feature of these sets of polynomials is that the reduction of any 
polynomial with respect to Q does not depend upon the order in which we use the polynomials 
therein in the reduction procedure. Another vital property is that given a monomial ordering, 
the Grobner basis (reduced with respect to itself) is unique for any given ideal. Unfortunately, 
one of the biggest hurdles in computational algebraic geometry is that the algorithm for 
determining the Grobner basis can be very intensive. 

B Primary Decomposition Algorithms 

In this section, we discuss in a little more detail the key algorithm used throughout the 
paper. There are now several primary decomposition routines available ^Hl 031201, many 
of which are implemented in algebra systems such as fOl ^] • We make extensive use of the 
algorithm due to Gianni, Trager, and Zacharias (GTZ) in this paper and so we shall 
now give a brief description of the basics of the algorithm's workings, following closely such 
texts as miEI]. 

The GTZ algorithm is built around the same splitting principle as was used in Subsec- 
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tionEm that is, if (/ : /°°) = (/ : /') for some /, then 



(40) 



Given this fact the GTZ algorithm works by specifying the polynomials / and by reduc- 
ing the primary decomposition of an ideal of dimension d to a problem involving primary 
decompositions of zero dimensional ideals. An existing algorithm can then be employed to 
primary decompose the zero dimensional ideals. We thus split our description into these two 
halves. First, we describe how the GTZ algorithm reduces everything to zero dimensional 
primary decompositions and finds a suitable /. Second, we give a brief discussion of how 
one obtains a primary decomposition of a zero dimensional ideal. 

B.l GTZ reduction 

The first step is to reduce the d dimensional decomposition problem to a dimensional 
one. We start with an ideal / in the ring C[Xi, . . . ,X„]. First, choose a maximal subset 
Y = {Yi, Y^} of the variables of the ring, X = {Xi, . . . , such that these variables are 
independent mod /. That is, /nC[Fi, Y^] = {0}. Geometrically, Y are the variables along 
L{I) and X\Y, transverse. Thus, d is the dimension of I. Now take the polynomials defining 
/ to be polynomials in Ic{y)[x\y] C C(Y)[X \ Y]. That is, pretend that the Y variables are 
coefficients. The ideal /c(y)[x\y], with all Y variables in / considered as coefficients, is then 
zero dimensional. 

Now, for our original ring C[X], choose a monomial ordering <, with Yi < Xj for all i 
whenever Xj G (X \ Y). Take a Grobner basis Q of I with respect to <. This is then also 
a Grobner basis of Ic{y)Ix\y], via restriction of < to X \ F. We are now in a position to 
isolate the / which GTZ employ. We take / to be the least common multiple of the leading 
coefficients of the poljTiomials in Q, with these polynomials taken to lie in C{Y) [X\Y]. The 
crucial observation is then the following: 



Thus, of the two halves of the saturation decomposition / = (/ : f°°) fl {I,f''), the first 
factor can be addressed by a zero-dimensional primary decomposition (to which we turn in 
the next subsection), leaving us with only /' = (/, /'), to deal with. We can then repeat the 
above process on I', and iterate until when there is nothing new in the second factor, i.e., 
when a factor we already have lies within the starting point for the next iteration. 



ic(Y)ix\Y] n C[X] = (/ : D . 



(41) 
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B.2 Zero dimensional Primary Decomposition 

Finding a primary decomposition of a zero dimensional ideal is relatively straightforward us- 
ing Grobner bases. Any zero dimensional ideal / can be put in a so-called "general position" 
with respect to the lexicographical ordering induced from Xi > ... > This is defined by 
the following properties: 

• The primes P{i) in the primary decomposition of / have a reduced Grobner basis with 
respect to the same ordering of the form 

{P{i)} = {Xi - /ii(X„), . . . , X„.i - /i„_i(X„), . (42) 

Here, we have hi G C[X„], i.e., they are simply polynomials in X„. 

• The ideals P{i) are coprime. In other words, the polynomials hi have as their greatest 
common divisors just an element of the coefficient field C, viz., a constant. 

Write Q for a corresponding minimal Grobner basis and define {h} = Q Ci C[X„]. There is 
then a theorem which states that ii h = h^^ ...h^j is the factorisation of h into a product 
of powers of pairwise non-associated irreducible factors, then the primary decomposition is 
just given by: 

/ 

An example of how this theorem can be used to implement an appropriate algorithm can be 
found in p/7j, as can various details. 

C Sturm Queries and Real Roots 

One of the topics of primary importance within this paper is the discussion of finding real 
roots of zero dimensional ideals. We shall thus briefiy describe some of the mathematical 
ideas involved in this appendix, following closely the excellent treatments of [211 122] • 

To commence, a finite set of polynomials within C[Xi, is zero dimensional iff any 

Grobner basis of the associated ideal contains a polynomial with leading monomial Xf' for 
each i G [1, fc]. Once a zero dimensional system has been identified one of the central notions 
in the study of its real roots is that of a Sturm query. Let P G M\x\ be a real polynomial 
and Z, a set of points. The Sturm query is given by the following expression: 

5Q(P, Z) = tl{x G Z\P{x) > 0} - tl{x G Z\P{x) < 0} . (44) 

If we had this function, then, for a zero-dimensional ideal / of real polynomials, the 
number of real roots is simply ^^(l, r(/)), where r(/) is the (discrete) set of real roots for /. 
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Moreover, we can also test sign conditions, another real algebro-geometric device which we 
use extensively in the paper. In such calculations we wish to know the sign taken by a given 
polynomial P evaluated at the elements of r(J). We note that, by definition, the following 
system of equations holds: 





1 


1 


1 ^ 


/ 







1 


-1 




v 





1 


1 ) 


V 



ti{xGr(/)|P>0)} 



/ 5Q(l,r(/)) \ 

5Q(P,r(/)) 
V SQ{P\riI)) ) 



(45) 



Once the Sturm queries are known, we can immediate solve for the quantities tl{x G r(/)|P = 
0, P > 0, or P < 0}, which are what we are after. One can also, in the same way, ask about 
the signs of lists of polynomials. This just involves the study of a bigger matrix equation. 

Thus we see that, once we know how to algorithmically compute Sturm queries, we can 
find the number of real roots of an ideal as well as the signs various polynomials take on 
those roots. How then is a Sturm query obtained algorithmically? The starting point here 
is to notice that if I is zero dimensional then the quotient ring Rq = R[Xi, . . . ,Xn]/I is a 
finite-dimensional P-vector space A. We can imagine taking a basis consisting of functions 
which are 1 on one root and zero on all the others, with one such function in the basis for 
each root. One can then obtain any function on the roots by combining multiples of these 
basis elements in the correct manner. We can define various linear maps on this space. One 
such map, Lf : A ^ A can just be defined to be multiplication within Rq by a function /. 
One can also consider bilinear maps Hg : A x A ^ R defined by Hg{fi, = Trace(L/j/2c,). 
Clearly the matrix associated to Hg in some basis for A is symmetric. 

A theorem due to Hermite states that the Sturm query SQ{g,r{I)) is simply given by 
the signature of this symmetric matrix. This is, in fact, intuitively obvious when thinking 
in terms of the basis described above. This matrix can be obtained algorithmically using 
Grobner bases [21] • Algorithmically the signature of symmetric matrices is easy to find. All of 
the eigenvalues of a symmetric matrix are real and are given by the roots of its characteristic 
polynomial. The number of positive roots is then determined by essentially Descartes' law 
of signs (or its generalisation, the Budan- Fourier theorem) [21], i.e., by examining the signs 
of the coefficients of the characteristic polynomial. 

The methods describe above are not necessarily the fastest way to obtain the results 
required, particularly the number of real roots [211 122] • They are however the simplest to 
understand. The reader interested in further details of these kinds of calculations is referred 
to \21\ 1^ . From a practical stand point, all of the algorithms concerned with real roots 
which we require have been implemented in jTT] by Tobis [22] • 
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